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We present a theory of the generation of magnetic seed fields in bubble collisions during a first- 
order electroweak phase transition (EWPT) possible for some choices of parameters in the minimal 
supersymmetric Standard Model. The theory extends earlier work and is formulated to assess the 
importance of surface dynamics in such collisions. We are led to linearized equations of motion with 
0(3) symmetry appropriate for examining collisions in which the Higgs field is relatively unperturbed 
from its mean value in the collision volume. Coherent evolution of the charged W fields within 
the bubbles is the main source of the em current for generating the seed fields, with fermions 
also contributing through the conductivity terms. We present numerical simulations within this 
formulation to quantify the role of the surface of the colliding bubbles, particularly the thickness 
of the surface, and to show how conclusions drawn from earlier work are modified. The main 
sensitivity arises such that the steeper the bubble surface the more enhanced the seed fields become. 
Consequently, the magnetic seed fields may be several times larger and smoother over the collision 
volume than found in earlier studies. Our work thus provides additional support to the supposition 
that magnetic fields produced during the EWPT in the early universe seed the galactic and extra- 
galactic magnetic fields observed today. 
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I. INTRODUCTION 

Identifying the source of the observed large-scale 
galactic and extra-galactic magnetic fields remains 
an unresolved problem of astrophysics One of 
the interesting possible sources is cosmological mag- 
netogenesis, where the seed fields would have arisen 
during one of the early- universe phase transitions. In 
this work our interest is seed field production dur- 
ing the electroweak phase transition (EWPT) dur- 
ing which the Higgs and the other particles acquired 
their masses. 

If magnetogenesis occurred during the EWPT, it 
most likely required a first-order phase transition. 
A first-order phase transition proceeds by a process 
in which bubbles of matter in the broken phase nu- 
cleate within the unbroken phase, similar to the fa- 
miliar process of steam condensing to water as the 
water- vapor mixture is cooled. Although it is gener- 
ally believed that there can be no first-order EWPT 
in the Standard Model d, there has been a great 
deal of activity in supersymmetric extensions [3], 



and for certain minimal extensions of the Standard 
Model there can be a first-order phase transition 0- 
6] . Limits on parameter-space of the minimal exten- 
sion of the standard model (MSSM) placed by elec- 
tric dipole moment measurements and dark matter 
searches allow a first-order EWPT which could lead 
to successful electroweak baryogenesis Q and the 
possibility that we are exploring, namely that seed 
fields responsible for the large-scale magnetic fields 
seen today are created during the era of the EWPT. 

Interest in these issues has led to quantitative 
studies of EWPT magnetogenesis based on the so- 
lution of equations of motion (EOM) derived from 
specific models. In the Abelian Higgs model, the 
first-order phase transition developed as the Uni- 
verse condensed into bubbles consisting of localized 
regions of space filled by Higgs field in a broken 
phase. This model was one of the earliest attempts 
to describe seed field production during the EWPT. 
The EOM related the seed fields to gradients in the 
phase of the Higgs field that were produced when 
bubbles merged following nucleation 

Simple 

and transparent solutions to the EOM evolved from 
specific field configurations applied at the point of 
collision in a relativistic 0(1, 2) symmetric model. 

The production mechanism of seed fields within 
an EOM approach has been pursued more recently 
within the framework of the MSSM [ill G3 along 
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the lines of the Abelian Higgs model in Ref. (8l-[To| 
in which 0(1,2) symmetric solutions evolved from 
specific field configurations applied at the point of 
collision. In this work, new EOM were derived from 
the MSSM Lagrangian, and accordingly the bubbles 
developed a coherent mode of charged W ± fields. 
As the bubbles merged, the charged gauge fields re- 
placed gradients in the phase of the Higgs field as 
the source of the electromagnetic em currents pro- 
ducing the magnetic seed fields, and the mechanism 
by which this occurred was developed in detail. Nu- 
merical results [12j showed that the MSSM produced 
seed fields of a size similar to the Abelian Higgs 
model even though the source of the current in the 
two approaches was quite different. 

Although these earlier studies gave insight into the 
production mechanism of seed fields in a first-order 
EWPT, applying boundary conditions at the time 
of the collision as implemented in the 0(1,2) for- 
mulations did not permit an assessment of the role 
of the dynamics of the bubble surface in the seed 
field generation. To explore the role of the surface 
it is necessary to specify the values of the W fields 
and their time derivatives on a surface (t = to, r, z) 
before the collision occurs. Thus, incorporating the 
initial stages of evolution of individual bubbles on 
the collision is an aspect of physics absent from the 
treatments found in Refs. |a-fl2|. 

In our more recent studies [H, [l4| results of calcu- 
lations were presented in which bubble surface dy- 
namics were taken into account. We identified a 
source of quantitative sensitivity to the bubble sur- 
face, and the results presented therein showed that 
the magnetic fields produced could be as large as, 
and possibly even larger than, those calculated in 
their absence. Encouraged by these results, in the 
present work we develop and extend the EOM for- 
mulation within the MSSM upon which this earlier 
work was based, along lines identified there. 

We begin, in Sect. II, by reviewing our EOM ap- 
proach. We also discuss briefly the nature of a first- 
order EWPT and some of the issues associated with 
the dynamics of the bubble surface that our theory 
is intended to address. In Sect. Ill we present our 
extended theory developed in 3 + 1 dimensions in 
a regime where the bubble collisions may be consid- 
ered "gentle" [HE2]. For gentle collisions the EOM 
linearize and display a transparent connection to the 
earlier work of Refs. pj-fiol [l2| . Various theoretical 
considerations necessary for assessing the role of the 
bubble surface in magnetic field generation, includ- 
ing the importance of establishing appropriate initial 
conditions, are developed. 



Out theory is applied in a specific model along 
the lines of Ref. [l2|, EH in Sect. IIV1 Because in our 
present formulation boundary conditions are applied 
before the collision occurs, we are able to examine in 
addition to collisions, the nucleation process, which 
is the evolution of the bubbles before the collision 
takes place. The numerical solutions of nucleation 
and collisions in this model are then presented in the 
following two sections. 

In Sect. |V] we examine the W fields, the current 
produced in collisions, and the magnetic field. For 
collisions, find the magnetic field to be larger in both 
scale and magnitude compared to our earlier 0(1, 2) 
results [l2j|. In Sect. IVll we estimate the sensitivity 
of the seed fields to the steepness of the surface of 
the scalar field, and in Sect. lVIll the sensitivity to the 
bubble wall speed and conductivity of the medium. 
Compared to results given in Refs. [3], we find that 
the magnetic seed fields are not only larger in mag- 
nitude but extend over substantially larger spatial 
scales than the results shown there. 



II. MAGNETIC FIELD CREATION 
DURING A FIRST-ORDER EWPT IN THE 
MSSM 

The EOM of this work are based on the same un- 
derlying MSSM Lagrangian as that of Ref. [HI . This 
Lagrangian is assumed to support a first-order phase 
transition and is of the form 

C MSSM = C l +£ 2 

+ lcptonic, quark, and supersymmetric 
partner interactions 

C 1 = -\wl v W^ v --B^B^ 

C 2 = \^d,- 9 -T-W,- 9 -B,W 

- V($,T) , (I) 
where T is the temperature and 

K> = W-d v Wl-ge ijk WlWi 

Bfj, u = d^B u — d u B^ . (2) 

Here W\ with i = (1,2), are the W+,W~ fields, 
$ is the Higgs field, and r l is the SU(2) generator. 
Fermions are not explicitly considered since earlier 
work to which we want to compare likewise ignored 
them. Because we are working within the framework 
of the MSSM, the bubbles that form as the phase 
transition progresses naturally consist of a region of 
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space filed by the Higgs field along with a cloud of 
other constituents of the MSSM Lagrangian in the 
broken phase. 
As in Ref. 



12j, we first derive "exact" EOM using 



an effective Lagrangian at the classical level from 
which the supersymmetric partners have been pro- 
jected out as explicit degrees of freedom, but whose 
effect is retained by a renormalization of the effective 
potential to maintain the properties of the first-order 
phase transition. These EOM are complicated non- 
linear partial differential equations coupling the W, 
B, and $ fields. From their solution one may obtain 
the physical Z and A em fields, 



A 1 



i 



Va 2 + a 
l 



- : (g'W^+gB, 



V 9 + 3 

In the picture we are developing, the Higgs field 
plays a central dynamical role in EW bubble nucle- 
ation and collisions, and therefore the effective (now 
appropriately renormalized) Higgs potential V($>, T) 
is an essential element in the theory. Although this 
potential is not known at the present time, depend- 
ing as it does on the unknown parameters of the 
MSSM as well as the properties of the plasma in the 
early Universe at the time of the EWPT, its spe- 
cific form is not relevant for the purposes of this 
paper. We require only that it should produce a 
first-order phase transition, consistent with certain 
MSSM extensions including for example those with 
a light right-handed Stop [15j . 

The various parameters are discussed in many 
publications 0] . For our calculations we use the lab- 
oratory values, 
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mw 



e/sm6 w = 0.646 
gta,n8w = 0.343 , 
80.4 GeV , 



(4) 



where mw is the mass of the bosons, and we 
define 



G = gg'/Vg 2 + 9' 2 = 0.303 



(5) 



In this section and throughout the paper units are 
such that H = c = 1, with distance and time ex- 
pressed in units of mw- 



A. First-order electroweak phase transition 

Coleman's model [l6| provides a conceptual 
framework based on a Lagrangian for understand- 



ing the phase transition. Although oversimplified in 
that it lacks medium effects, which can lead to an 
asymptotic wall speed v wa u < 1 depending on the 
pressure difference in the true and false vacuum, it 
is seminal in that it was one of the earliest EOM 
descriptions of the physics of a first order phase 
transition. In his model, prior to nucleation, the 
dynamics of bubbles is formulated in the Euclidean 
metric, which is 0(A) symmetric. After nucleation, 
the bubble expands in 0(1, 3)-symmetric Minkowski 
space-time [16| . As the phase transition develops the 
bubbles start to merge or " collide" . Eventually they 
completely merge, at which point the phase tran- 
sition is completed. In subsequent work magnetic 
fields have been understood to be created in the col- 
lisions of bubbles. 

In the Lagrangian of Eq. (fT]), the EWPT is driven 
by nucleation of the scalar Higgs field just as imag- 
ined in Coleman's model [l6j]. In this picture, the 
vacuum state of the Universe corresponds to a local 
minimum in V($,T). The phase transition occurs 
as the temperature T is lowered through the tran- 
sition temperature T c when V($,T) develops a de- 
generate second minimum at a larger value of < $ > 
separated from this minimum by a barrier. As the 
universe continues to expand and cool, the depth 
of the second minimum increases, meaning that the 
Universe can lower its energy by moving from the 
original, now metastable, false vacuum to the lower 
energy true vacuum. Because the two minima are 
separated by a barrier, the transition from the false 
to the true vacuum is delayed as the temperature 
continues to drop, a process referred to as super- 
cooling. This delay influences bubble characteris- 
tics, and a first-order phase transition is accordingly 
classified as weak or strong depending on the degree 
of supercooling. A comprehensive phenomenologi- 
cal study of the kinetics of cosmological first-order 
phase transitions, such as the EWPT, in terms of 
such an effective potential is given for example in 
Ref. (H. 



B. Bubble dynamics and magnetic field 
creation with 0(1, 2) Symmetry 

In their analysis^ of the Abelian Higgs model, Kib- 
ble and Vilenkin [8[ obtained magnetic fields as bub- 
bles merge in a regime of gentle collisions. They 
obtained EOM in this case by making an expan- 
sion about point p(x) = po (the "Kibble- Vilenkin 
point"). From these EOM, expressed in terms of 

), where t is the time 



the variables (z, r = \[W 
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and r — \/ x 1 + y 2 is the distance of a point from the 
z-axis, they obtained 0(1, 2) symmetric solutions us- 
ing jump boundary conditions applied at the time of 
collision. They demonstrated that when the phase 
of the Higgs fields is initially different within each 
bubble an axial magnetic field forms as the bubbles 
merge and that this field has the structure of an 
expanding ring encircling the overlap region of the 
colliding bubbles. 

In our earlier work in the MSSM [l2|, EOM were 
obtained by making an expansion about the Kibble- 
Vilenkin point. First, p(x) was expressed as 

p(x) = Po + aSp(x) , (6) 

with po the magnitude of the mean scalar field at 
the center of a single bubble and aSp fluctuations 
of the magnitude in the scalar field once the bubbles 
merged. Making an expansion in a as in the Abelian 
Higgs model we obtained linearized equations within 
the bubble overlap region. Collisions in which the 
Higgs field is relatively unperturbed from its mean 
value when the bubbles merge were termed "gentle." 

Then, assuming as in the Abelian Higgs model 
that the collision begins at time t = to (called t c 
in Rf. [Hj]): when the bubbles first touch at z = 
0, we used jump boundary conditions to determine 
the charged W fields in 0(1,2) symmetry and the 
magnetic field. These boundary conditions recognize 
that in some collisions the sign of the z-components 
of the W + field at leading order in a is opposite 
in the two colliding bubbles while the z-components 
for W~ has the same sign, and that for others the 
phases are the same for in the two bubbles. For 
the first case, referred to with a superscript I, the 
boundary condition was 

w zI (T~to,z) — w e(z) 
^-w zI {r = t ,z) = 0, (7) 

where e(z) is the sign of z, and for the second, iden- 
tified with a superscript II, 

w zII (t = t ,z) = w 
^-w* II (T = t ,z) = . (8) 

Comparing to the Abelian Higgs model we found the 
two magnetic fields to be of similar size. 



C. Effects of the medium: bubble surface 
motion and conductivity 

Because effects associated with the surface break 
0(1,2) symmetry, a theory formulated with this 
symmetry and the associated (z,t) variables such 
as that of Ref. 12] is unsuitable for exploring the 
consequences of surface dynamics. To explore such 
effects, the theory has to be formulated in 3 + 1 di- 
mensions, and the appropriate symmetry group is 
0(3). 

Additional drawbacks of 0(1, 2) symmetry include 
the restriction that v wa ii — c and difficulty includ- 
ing electrical conductivity in Maxwell's equations. 
Values of v wa u < 1 have been obtained in diverse 
studies including modeling collisions of bubbles with 
constituents of the plasma [l8| and solving EOM 
including nonlinear terms based on an MSSM La- 
grangian [l9j . 

A value v wa u ^ 1 and finite conductivity both 
directly affect the seed field. This has been dis- 
cussed comprehensively and their effect estimated 
in Refs. [a-tlPl] within the context of the Abelian 
Higgs model. It was found there that finite con- 
ductivity would lead to the decay of the currents 
(and therefore the magnetic field) with a character- 
istic time td ~ a/m [8j with m the gauge boson 
mass. An additional consequence of the large con- 
ductivity arises as follows. Since the magnetic fields 
propagate with the speed of light, for slowly expand- 
ing bubbles these fields would very quickly escape 
from the region of the bubble collision in the ab- 
sence of conductivity. However, because of the large 
conductivity the magnetic fields become "frozen" or 
confined to the region of the bubbles, hindering the 
escape of magnetic flux into the surrounding false 
vacuum. Kibble and Vilenkin showed that the loss 
of flux is negligible provided that aR c v >> 1, where 
R c is the bubble radius at collision time. With val- 
ues of conductivity that are believed to characterize 
the plasma, currents and magnetic fields persist on 
time scales that are long compared to those of the 
symmetry breaking scale. 



III. EQUATIONS OF MOTION IN THE 
MSSM WITH 0(3) SYMMETRY 

In this section we develop a general framework 
in O (3) symmetry extending our earlier work in 
Refs. [ll GJ] based on the MSSM. Our formulation 
is intended to be capable of following the evolution 
of bubbles with a given wall speed v wa ii in 3 + 1 
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dimensions, starting at time of nucleation, and de- 
termining the magnetic field generated in collisions 
including effects of finite conductivity. 

We begin the development of our theory in 
Sect. IIII Al by extending the concept of a gentle colli- 
sion to the case where the surface is explicitly consid- 
ered. We then derive, in Sect. IIII Bl EOM by making 
an expansion of the scalar field for a pair of bubbles 
about the mean scalar field. The expansion, justified 
for gentle collisions, leads to linearized EOM, thus 
simplifying the theory. In Sect IIII Cl we discuss some 
of the new issues that are encountered in solving the 
EOM when the bubbles are initially separated. As 
shown in Sect. IIII Dl the same expansion leads to an 
expression for the em current in terms of the W and 
to a corresponding Maxwell equation. 

A. Gentle collisions in electroweak theory 

When the surface is considered, we generalize 
Eq. ([5]) by writing 

p(x) = p(x) + aSp(x), (9) 

with p{x) a simple function approximating the mean 
scalar field at any point x in the medium in the col- 
lision. The quantity aSp(x) is, as above, the change 
of the magnitude in the mean scalar field induced 
by the collision. In this paper we will obtain, as in 
Ref. [l2|], linear approximations to the exact non- 
linear EOM by expanding them in terms of the pa- 
rameter a appearing in Eq. ©. The resulting EOM 
are similar to those of Ref. [12| , but because we now 
have the surface to consider they differ in a num- 
ber of essential ways and require the development of 
completely new techniques to solve. 

The justification of the expansion in terms of a in 
the present case arises as follows. Clearly, when two 
colliding bubbles are completely separated p(x) for 
these bubbles is, to a very good approximation, the 
sum of the scalar fields of independent bubbles, and 
there are no significant fluctuations that need to be 
considered (a ss 0) to the extent that one has con- 
fidence in the choice made for the scalar field in an 
individual bubble. Additionally, for two completely 
interpenetrating bubbles p(x) within the central re- 
gion is approximately the same as the scalar field at 
the center of a single one of the colliding bubbles, 
and the justification of the expansion is the same 
there as it was in Ref. 

The new issue is to justify the expansion in the 
peripheral region when two bubbles first begin to 
merge. The critical point to recognize is that in this 



region the size of a and hence the accuracy of the 
expansion will depend on how well the mean field 
p(x) that appears in the EOM approximates the ex- 
act scalar field for the colliding bubbles. We come 
back to this issue below. 



B. EOM in electroweak theory for gentle 
collisions 

The fact that i/j and W d (for d = (1,2)) enter 
quadratically in the p-equation (Eq. (8) of Ref. [l2j ) 
places two important constraints on these quantities: 
(1) ip and W d must have an expansion in odd powers 
of a 1 ' 2 , if we require the square of these quantities 
be analytic in a; and, (2) expanding this equation 
to leading order in a 1 / 2 , we find that the terms V > 
k/ ) 1 , and u/°) 2 must vanish. This is most easily 
seen in the Euclidean metric, from the fact that the 
square of each enters with the same sign. However, 
the same must be true in the Minkowski metric as 
well by analytic continuation. In view of these con- 
siderations, ip v and Wj~ for d = (1,2) have the fol- 
lowing expansion 

^{x)=a l ' 2 ^+a 3 ' 2 ^ + ... (10) 

W d = a^ 2 w^ d + a 3 ' 2 w^ d + ... . (11) 

It is natural that an expansion in the same param- 
eter a 1 / 2 remains appropriate for d = 3. However, 
there is no requirement that the leading term vanish, 
so we take 

W 3 = w^ 3 + a l ' 2 w^ 3 + a 3 / 2 w^ 3 + ... . (12) 

With a « we may take p(x) ~ p(x), and the 
B-, 0-, and W^-equations then give, to first order in 
a 1 / 2 , 

[d 2 + ^f(g 2 + 9 ' 2 )}^=0 (13) 
d a ^ = (14) 

where now 

^(x) = d a Q- {9 +9 2 }/ Zg). (15) 

Equations for w^ d may be obtained by expanding 
the B- and W^-equations through order a 1 / 2 . For d = 
1 or 2 (corresponding to d' = 2 or 1, respectively), 
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we obtain the pair of equations 

= d 2 w u ^ d - d v d ■ + m{xfw^ d 

- 2[d^(w^ 3 w^ d ' -w^ d 'w^ 3 ) 

- (wW d 'd u w^ 3 - w^d v w^ d ')\ 

- 4[(w^ 3 ) 2 w^ d - w^ 3 ■ w^ d w^ 3 } , (16) 

where m(x) is the mass of the W field, 
p{xfg 2 



m(xY 

The corresponding equation determining w. 
d = 3 is 



(17) 



(i)d 



for 



d 2 w £ )3 - a u a ■ w {1)3 = p(x) 2 gi>P , (is) 

which can be solved once the driving term if)^'(x) 
has been independently determined from the solu- 
tion of Eqs. (|13I14|) . The considerations for fixing 

the boundary conditions for w^ d and wi^ 3 are sim- 
ilar and discussed below. 

This field w^ 3 appearing in Eq. (1161) is found to 
be the solution of 



^(0)3 _ . w (0)3 







(19) 



Because no mass appears in this equation, W occu- 
pying this mode propagate at the speed of light and 
experience no interaction with the scalar field of the 
bubble to lowest order in a, unlike the W described 
by w v . Because of this, there is no appreciable 
coupling to wj? 3 , and Eq. (TT6"]) becomes 



= d 2 w a v - d v d ■ w a + m{xYw a v . 



(20) 



We see that for sufficiently gentle collisions, all rele- 
vant equations are linear in W . 

Simplifying the non-linear p equation (Eq. (8) of 
Ref. [ijj) using the fact that w" and ip v in leading 
order go as a 1 / 2 , we find that to leading order in a 
p{x) satisfies the equation 







d 2 p(x) + p(x) 



dV 
dp* 



(21) 



and the solution of this equation is clearly identified 
with p[x) appearing in Eq.®. Methods for solving 
Eq. (|21l) are discussed in many places, for example 
Ref. [23. 

In 0(3) the complete set of EOM for describing 
bubble collisions with the surface considered has now 



been derived and consists of coupled partial differen- 
tial equations for the relevant fields. Strictly speak- 
ing, the set is non-linear because the solution of 
Eq. (|2"Tj) is coupled to the W field through the mass 
of the W, evident in Eq. (TTT)) . 

However, because Eq. (|2T|) does not depend on 
, the coupling of the magnitude of the scalar field 
p(x) to the charged W fields w^(x) in Eqs. (|20I21|1 
is particularly simple and allows p(x) for a system 
of colliding bubbles to be determined once and for 
all. Equation (f2"U|) is then effectively uncoupled from 
Eq. ([2~Tj) and may be solved directly to obtain w% (x) 
for all x, effectively linearizing the EOM and result- 
ing in an enormous simplification. That the EOM 
are effectively linear implies that for gentile collisions 
the coupling of the W to the Higgs dominates the 
self-coupling of the W fields. 



C. Surface effects and W ± fields in bubbles 



With the solution of Eq. (|2Tj) nearly constant at 
p(x) — pa in the broken phase comprising the in- 
terior region of single or overlapping bubbles, our 
previous work [12| was simplified. However, now 
that we are considering as well the surface, where 
p(x) begins dropping to its value in the symmetric 
phase p{x) = outside we cannot ignore the spa- 
tial dependence of the mass as given in Eq. (fT7|) . 
The spatial dependence in the surface not only in- 
troduces a few technical challenges but also, as we 
will see, new physics with quantitative significance 
not present in [12j . 

One of the consequences of the spatial dependence 
of the W mass can be seen by taking the four- 
divergence of Eq. (|2"0")) . By so doing, we obtain the 
auxiliary condition 



X a (x) = 



(22) 



where 



cV(m 2 w a -dm 2 
X a {x)^ 1 2 ^ =d-w a + W ° . (23) 
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Equations (|22I23I) require 



d-w a = 



(24) 



Thus, in contrast to the calculation in Ref. [12j, it is 
no longer true that d ■ w a = 0, and as a consequence 
we find that the equations of motion for the W fields 
more complicated. 
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The physics becomes clearer by using the relation- 
ship in Eq. (J24J) to rewrite the EOM in Eq. ([20]) as 



= d 2 w" 



<j a ■ dm 2 



+ m 2 w° 



(25) 



The solution to this set of equations is equivalent to 
the set in Eq. (|2"0")) provided the auxiliary condition 
Eq. (|22|) is maintained for all (t,x). 

The transformed EOM Eq. reveal that the 
spatial dependence of the mass provides a 

perhaps unexpected sensitivity to the bubble sur- 
face. The sensitivity occurs through the term w a ■ 
dm 2 /m 2 , which becomes in fact divergent in the 
limit of an infinitely sharp bubble surface. At this 
point one cannot rule out significant modifications 
to results obtained in 0(1,2), where the surface is 
ignored. 

To see how the auxiliary condition Eq. (|2"2"|) may be 
maintained for all (t,x), note that Eq. (|25p requires 
X a {x) to satisfy the Klein-Gordon equation 

d 2 x a (x)+m(x) 2 X a (x) =0. (26) 

By choosing the initial configuration of w a {x), at 
time t — to, to satisfy 



x a {h,x) = o 



and 



dx a (to,x) 
dt 



, 



(27) 



(28) 



we assure that x a { x ) = f° r ai l future times since 
Eqs. (|27|28|) are boundary conditions for the trivial 
solution X a (t, x) = of Eq. Thus, Eqs. (12TI25)) 

provide constraints on the initial conditions for the 
W ± fields. 

To establish the initial conditions requires the 
choice of a time to at which the initial values of the 
W fields in the bubble are specified. In Ref. [12| the 
counterpart of to was the point of first contact of 
the bubbles. In the current approach to may be in 
fact much earlier, in particular it could be as early 
as the time of nucleation t n . The choice of initial 
conditions is further discussed in the context of a 
our model in Sect. IIV E 31 below. 

These observations make it natural to distinguish 
two categories of initial conditions when the surface 
effects are considered. The first, which we will refer 
to as boundary conditions, consists of the initial W 
fields that may be chosen freely. The second consists 
of the set determined by Eq. (|27I28[) . which we will 
refer to as the constrained initial conditions. 

The definition of x a given in Eq. (|23p requires 
m(x) 2 > everywhere as it would be in the mean- 
field approximation adopted in Sect. IIVDI 



D. Maxwell's equations 

We may find the Maxwell equation for the em field 
Af, m (x) by taking the linear combination of the 
and B indicated in Eq ©. An expression for the 
corresponding em current j v (x) consisting of terms 
quadratic and cubic in the three fields W l {x) imme- 
diately follows !12|. 

The result for j u {x) may also be simplified by ex- 
panding the A em and W fields in powers of a 1 / 2 . 
Letting (x) refer to the terms in the expansion 
of Al m (x), we find that to leading order Af, m (x) — 

(2) 

al ' (x) and satisfies the following Maxwell equation, 

d 2 a^ - d v d'aF> 

= Ge ab3 (w^ b d-w^ a 

- w^ a d v w^ b + 2w^ a ■ dwP b ) 

= iTrjrix) , (29) 

where we have introduced the coupling parameter G 
defined in Eq. ([5]). From this we learn that the first 
non-vanishing contribution to the em current is of 
order a 3 / 2 and that it depends on the components 
w^ 1 of the charged W fields (i = 1 and 2), calculated 
at order a 1 / 2 . Expressing the current in terms of G, 
we find 

4TTj* m (x) = Ge ab3 (wl 1)b d-w {1)a 

- w^ a d I/ w^ b + 2w {1)a ■ dwi 1)b ){30) 



It is easy to prove that this current is conserved, 



d-j em (x) = 



(31) 



using the fact that at the classical level 
[w a {x),w b v {x')] = and the fact that the W 
fields appearing in Eq. (|3"0")l satisfy the EOM, 
Eq. O. 

So far the bubble has been considered to consist 
purely of the scalar Higgs field and the associated 
cloud of charged W gauge bosons coupled to it. Ob- 
taining the contribution of the charged gauge 
fields, Eq. (|30|) . to the em current for collisions of 
such bubbles is one of the important results of this 
derivation. 

However, fermions also contribute to the em cur- 
rent and have a significant impact on magnetic seed 
field production. One contribution was discussed re- 
cently in Ref. [l9| and estimated there for the nucle- 
ation phase of the collision. Another, discussed in 
Sect. Ill CI occurs through the conductivity of the 
medium a. This is one of the most important and 
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best-known contributions and may be taken into ac- 
count through its associated current j c (x), 

j c (x) = aE(x) , (32) 

where the usual assumption that j c (x) is propor- 
tional to the electric field E has been made. De- 
tailed calculations of a in the early universe are avail- 
able [2l|. 

To find Maxwell's equation for the magnetic field 

B, 

B = V x A em , (33) 

we multiply Eq. (|29l) by eij k dj, obtaining 

e l3k d 3 d 2 Al m - e ijk djd k d ■ A em 

= e t3k d 3J r . (34) 

Expresing Eq. in components, 

Bi = e l]k d 3 A e k m , (35) 
we immediately find the desired result, 

8 2 B = 5 x j em . (36) 

IV. MODELING BUBBLE COLLISIONS IN 
0(3) SYMMETRY 

We will assess the importance of surface effects 
by making numerical simulations that can be mean- 
ingfully compared to earlier work in Refs. [HI, [l4j]. 
The common dynamical framework is summarized in 
Sect. IIV A[ with common geometrical aspects spec- 
ified in Sect. IIVBI The representation of the mean 
scalar field for two bubbles p{x) — pi(x) in this ge- 
ometry, along with the arguments for its choice, is 
given in Sect. IIV CI The presence of other bubbles 
are taken into account in a mean-field approximation 
in llV Dl Initial conditions on the W fields, including 
boundary conditions and the constraints imposed by 
surface geometry, are discussed in Sect. IIV El We 
specialize the theory to cylindrical coordinates in 
Sect. IIV Fl The model is applicable to both weak 
and strong first-order phase transitions and incor- 
porates some important medium effects absent from 
our former work. 



A. The dynamical framework 

The familiar conceptual features of our MSSM 
theory have been presented in Sect. Ill Al When the 



Higgs couples to the other fields as it does in the 
MSSM through the Lagrangian of Eq. ([T]), strong 
and highly non-perturbative couplings arise forming 
a tightly coupled many-body system as the phase 
transition develops. Specific consequences of this 
were suggested and embodied in Ref. [12| , and these 
apply as well to the present work. 

One of the identified consequences of this coupling 
is that as the bubble growth occurs, the (and 
other constituents that we are ignoring here) that 
enter the bubble from the plasma gain their mass at 
the expense of thermal energy, a cooling process that 
continues as volume available to the W increases 
with the bubble expansion. Another is that the cou- 
pled fields tend to follow the evolution of the Higgs 
field coherently. As they lose thermal energy, the W 
passing into the bubble enter a single mode, a solu- 
tion to the EOM discussed in the previous section. 

This mode plays a special role, and it is quite 
different from the more familiar incoherent thermal 
modes outside the bubbles. One may think of it as 
a coherent state, much like a state of electrons in 
a superconductor (except that the W are bosons). 
The mode of course evolves in time according to the 
EOM, and as the bubble expands and displaces more 
of the volume of the plasma the occupation of this 
mode also grows. The dynamics driving it is clearly a 
non-equilibrium component of the phase transition, 
and it is a basic assumption of our EOM approach 
that these coherent fields give rise to the seed fields 
as bubbles merge before thermal equilibrium is re- 
established. 

There are of course many such field configurations 
that satisfy the EOM, and the one that is realized 
in a given bubble in the phase transition depends on 
the overall history of the process just discussed. To 
calculate the net magnetic field produced in bubble 
collisions properly, one would have to evaluate the 
field corresponding to each possible initial configura- 
tion and average over the ensemble of configurations. 

The net effect of this averaging procedure was ex- 
amined in Ref. |12j and found to be factor less than, 
but the order of, unity. Thus, to get a fair estimate 
of the net magnetic field it is not necessary to explore 
the full range of possible initial conditions; rather, 
it is sufficient to examine one initial condition char- 
acteristic of the entire ensemble of possibilities. 



B. Bubble collision geometry 

For the application of our theory to the collision 
of two bubbles, we will be assuming, as in Ref. (l4| . 
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that the bubbles are nucleated simultaneously at 
time t — t n at points z — ±zq located symmetrically 
about the origin on the z-axis with nucleation radii 
r ns . Additionally, here we assume that the bubbles 
are well separated and non-overlapping before they 
collide. 

The scalar field of a single bubble may be repre- 
sented in general as 

p(t,x) = P J s (t,x) , (37) 

where f 3 (t, x) is the shape of the field and p c is its 
magnitude at the center. The f s (t,x) that we will 
use in this paper, essentially equivalent to that in 
Ref. [HI, is defined in Appendix [S] in terms the bub- 
ble nucleation time t n , nucleation radius r ns , surface 
speed v wa ii, and surface diffuseness a s . The surface 
diffuseness is approximately half the distance over 
which the scalar field falls from its 10% to 90% at 
t n for large bubbles, i.e. bubbles with r ns > 2a s . 

The bubble nucleated at z = z > is referred to 
as the right-hand (R) bubble and the one nucleated 
at z = — z < as the left-hand (L) bubble. We 
assume axial symmetry so that the relevant spatial 
coordinates (r, z) for a point are its axial coordinate 
r = \J x 2 + y 2 , its distance from the z-axis, and its 
longitudinal coordinate z, its distance from y — z 
plane that passes through the origin at z — 0. 

It is helpful to think of the collision in terms of the 
evolution of spatial surfaces that separate regions 
of the collision occupying true vacuum from those 
of false vacuum. The connectedness of the surfaces 
change as the collision process proceeds. 

Initially, for times t n < t < t c , where t c is the 
collision time or time of merging of the bubbles, the 
boundary consists of two disconnected surfaces, one 
for the left-hand bubble i — L and the other for 
the right-hand bubble i = R. For times t > t c the 
bubbles coalesce to form a region i = c with a single 
boundary surface, S c . 

The radius of the bubble surface R\/ 2 (t) is defined 
as the distance from the center of the bubble (its 
nucleation point) to the point at which the scalar 
field has fallen to half its central value. It of course 
depends on our choice of scalar field, whose details 
are given in Appendix [X] The collision time may 
then be taken to be the time at which the radius of 
either bubble becomes equal to the distance of its 
nucleation point from the origin of the coordinate 
system, i.e. when R\/2{t c ) = zq. The solution of 
this equation for our scalar field is given in Eq. (IA9|) . 

The spatial points forming the boundary surfaces 
are then determined by the equation 

R(r,z) = R 1/2 (t) , (38) 



where 

R(r,z) ee R R (r, z)6{z) + R L (r, z)6{-z) , (39) 

with 

R R (r, z) = Vr 2 + (z - z Q ) 2 (40) 

being the distance from the center of the right-hand 
bubble at zq to the point (r, z) on the surface, and 

R L {r, z) = ^r 2 + (z + z Q ) 2 (41) 

being the distance to the same point from the center 
of the left-hand bubble at — zq. 

After the bubbles merge, the right-hand and left- 
hand surfaces intersect on the x — y plane in a circle 

of radius bit) = J Ri/2{t) 2 — z q that expands at a 
rate 

db(t) = R 1/2 (t) dR 1/2 (t) 
dt b(t) dt [ ' 

Our model of scalar field in the collision of two bub- 
bles, including the region of coalescence is given in 
Sect. HVCl 

Since the bubble collision geometry is axially sym- 
metric with the relevant spatial coordinates being 
(r, z), where z is the distance of a point from the 
y — z plane on the z-axis, and r = \/ x 2 + y 2 being 
its distance from the z-axis, cylindrical coordinates 
is the natural coordinate system for expressing re- 
sults. In this coordinate system, the W fields may 
be taken to have the following form, 

wl{t,r,z) = w%(x), v = Q 
w a v {t,r,z) = x u w a {x), i/ = (l,2) 
v£(t,?,z) = w a z (x), u = 3, (43) 

with m° = — w za and Wq — w 0a . Correspondingly, 
we write the em current in cylindrical coordinates as 

jv(t,r,z) = (j (t,r,z),rj(t,r,z),j z (t,r,z)) . (44) 

C. Mean scalar field in two colliding bubbles 

An expression for the mean scalar field p(x) = 
P2{x) for a pair of colliding bubbles may be expressed 
in terms of the non-overlapping portion of the scalar 
field, Ap(x), 

Ap(x) = p L (x)(l - p R (x)/p c ) 

+p R {x)(l -p L (x)/ Pc ) , (45) 
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where the scalar fields in the left-hand and right- 
hand bubbles, while still separated, are taken as 

p L (t,r,z) = p c f s (t,r,z + Zq) 
p R (t,r,z) = p c f s (t,r,z - z Q ) (46) 

with f s (t,r,z) and p c given in Eq. (|37|) . 

Now, if the scalar fields were to simply add in the 
overlap region (which we know is not the case [Hj]) 
the scalar field would be 

Pl(x)+pr(x), (47) 

which is as large as twice the size of that when they 
are well separated. What we have learned from the 
calculation in Ref [12[ is that the scalar field, when 
the bubbles have completely merged, is actually half 
this value, so if we define 

p 2 (x) = Ap(x) + -(p L (x) + PR {x) - Ap{x)) , (48) 

then we obtain our ansatz, which is in agreement 
with the results of Ref. [14j . 

Because there is not much a priori guidance on 
how to define p{x) in the region where the bubbles 
begin to overlap, the accuracy of the expansion in 
terms of a appearing in Eq. @ is accordingly dif- 
ficult to establish for any given p(x). However, be- 
cause the peripheral region extends over a small vol- 
ume relative to that of the two bubbles, it is likely 
that the corrections for any reasonable choice will be 
relatively insignificant. 

The evolution of the scalar field in a collision with 
pi (x) defined in Eq. ([4*5)) is shown in Fig. [TJ This fig- 
ure shows the scalar field on which the calculations 
of Sect. [V] are based. The individual bubbles nucle- 
ate at time t n — on the z-axis centered at z — ±z , 
with zq = 35 and nucleation radius r ns = 20. Their 
scalar fields taken to have wall thickness a s = 4 and 
Vwaii = 1- Aside from our choice of surface diffuse- 
ness and some compensating adjustment in r ns and 
zq (the parameters zq and r ns are a bit larger than 
the ones used there so that with the thicker surface 
the bubbles do not overlap at t = to) the parame- 
ters and collision geometry closely matches that of 
Ref. [lif . Because of the larger surface diffuseness, 
this corresponds to a somewhat weaker phase tran- 
sition than that of Ref. [14j and thus provides an in- 
teresting contrast to the calculation presented there. 

The collision time is found from Eq. (|A9|) to be 
t = t c « 11. The bubble radii at this time are 
R\/2{tc) = zq = 35. For times t < t c the scalar field 
is approximately confined to two isolated regions, 
the two individual bubbles, and that after this time 
to just one region, the collision region. 




FIG. 1: Evolution of scalar field in the collision of two 
bubbles. Nucleation is at t n = 0. Field is plotted as a 
function of z for r = over the time interval < t < 40. 



D. Medium containing many bubbles 

In a medium undergoing a first-order phase tran- 
sition there are many distinct regions Si of scalar 
potential, each of which is either a single bubble or 
a collection of mutually overlapping bubbles. The 
net scalar field p(x) for this system may thus be ex- 
pressed as a sum over the scalar fields p^ l \x) for 
each of the N regions, 

JV 

p(x)=Y^P (i) (x) ■ (49) 

i 

For such a system, the mass of a W boson at any 
point x continues to be given by Eq. (fTT)) . where p(x) 
is the net scalar potential at the location of the W . 

In this paper we are interested in the evolution of 
at most two bubbles in this sum. But in tracking 
their evolution, we should not ignore the presence 
of the other bubbles. We will account for them in 
a mean-field approximation by averaging over their 
locations and representing their collective effect by 
an average scalar field p avi 

N' 

i 

N' 

= yJd^X^Ht.x) , (50) 

i 

where V is the total volume over which the integral 
runs and the prime on the sum means we exclude 
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the region(s) of explicit interest in the average. The 
quantity p av thus acquires its specific value from the 
presence of other bubbles that appear as the phase 
transition develops. 

Since p av is essentially the average scalar field aris- 
ing from the bubbles in the medium at the time of 
interest, at the onset of the phase transition p av rs 
and at the completion of the phase transition the 
scalar field becomes uniform with p av — > pq. Thus, 
Pav I ' Pa, with < Pav I Pa < 1, is not only a mea- 
sure of the extent to which the phase transition has 
evolved but also tracks the relative density of the 
bubbles in this evolution. One might expect that 
the binary collisions that are of most interest in this 
paper dominate for p av /po ^ 0.1 — 0.2 and that for 
larger values simultaneous collisions of multiple bub- 
bles begin to contribute significantly. 

Accordingly, taking into account the average value 
of the mean scalar field p av , and as long as p av /po ^ 
0.1 — 0.2, the variation of the mass of a W in the 
left-hand and right-hand bubbles when they are sep- 
arated is given by 

g 2 

m L {x) 2 = —{pav+PL{x)) 2 

g 2 

m R (x) 2 = — (p av +p R (x)) 2 , (51) 

respectively. In general, for a pair of bubbles in a 
medium experiencing any degree of overlap, m(x) is 
given by 

m{xf = 9 -{p2{x) + Pav f , (52) 

where pi{x) is defined in Eq. (|48l) . 

With m(0) referring to the mass at the center of 
one of the regions Si of the collision, clearly 

™(0) 2 = = m w ■ ( 53 ) 

Since it is also true from Eq. (|52"|) that 

™(0) 2 = y(p 2 (0)+p a ,) 2 , (54) 

and the scalar potential pi (0) = p c it follows that at 
the center of a bubble in the medium the scalar field 
is 

Pc = Pa- Pav ■ (55) 

The quantity p av determines the mass of the W in 
the bubble surface and thus plays a role in confining 
the W to the region of the bubble. Beyond this, the 
results of the theory are relatively insensitive to the 
specific value of p av - 



E. Initial conditions on W fields 

A meaningful comparison of the present theory 
to those of Refs. [HI, [l4| clearly requires that similar 
initial conditions be applied. In contrast to Ref. [Hj], 
the initial conditions are imposed here on W fields 
in bubbles that are initially separated. 

Although our intent is to assess the importance 
of surface dynamics by comparing the results of our 
present paper to those of Refs. [H, [Hj], we should 
keep in mind that there is no well-defined limit in 
which the present theory, which exhibits Q(3) sym- 
metry, exactly coincides with that of Ref. [12J , which 
exhibits 0(1,2) symmetry. 

1. Boundary conditions for individual bubbles 

As indicated earlier, in the dynamical framework 
of Refs. [HI, [l4| as a bubble expands it displaces 
W in the thermal plasma which enter the bubble in 
a single mode. The normalization of this mode in 
the bubble (that is, in the true vacuum) is deter- 
mined by requiring that the average number density 
of in the bubble be equal to the number density 
of the W quanta in the displaced plasma. This 
condition maintains the average density of W + to 
be roughly constant as a function of time and equal 
to the density of the W~ inside the isolated bubble. 
As a consequence 

\w^ 1 (t,x)\ = \w a = 2 (t,x)\ (56) 

with random relative phases. The linearity of the 
EOM in Eq. (l2Cfl) guarantees Eq. (|56| in isolated 
bubbles if a some initial time to w° has the same 
magnitude and shape for both a = (1,2), 

w^ito^x) = aw a = 2 {t ,x) 

J^«=i( t0) £) = a^-w a = 2 (t ,x) , (57) 

where a is an arbitrary phase. 

Because the electromagnetic current in Eq. (pOj) 
is antisymmetric in the labels a and b this current 
will vanish since the field w a (x) for a = 1 has the 
same dependence on x as that for a = 2. Thus, in 
an isolated bubble the electromagnetic current will 
vanish. 

However when two bubbles, each meeting the con- 
dition in Eq. (|56|) collide, it is in general not the 
case that w^{x) satisfying the equations of motion 
will be proportional in the region of overlap, and as 
a consequence em currents will form. This is the 
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reason why magnetic fields are in general produced 
when bubbles collide. We model this below by solv- 
ing the equations for the W helds with different sets 
of boundary conditions on w^, for a — 1 and a = 2. 



2. Boundary conditions for bubble collisions 



With Eq. (1571) in mind, we turn our attention to 
finding boundary conditions for collisions so that 
the W ± fields are as similar as possible to those of 
Refs. PUTS. By 

requiring that the bubbles be sepa- 
rated at the initial time t = to [HI and that Eq. (|57)l 
is satisfied there is no field produced until the col- 
lision occurs, as in Refs. 

The analog of the jump boundary conditions in 
Eqs. (EGO) at t in 0(3) are then 



and 



wl(t ,r,z) 
dwijto,^) 
dt 

dwf J fto,x) 
dt 



w L (r,z) - w R (r,z) 
dw L (r,z) dw R (r,z) 



dt 



<)1 



(58) 



w L {r,z) + w R (r,z) 
dw L (r,z) dw R (r,z) 



dt 



(59) 



The functions wl{x) and w R (x) are clearly the pro- 
files of the W fields in the left-hand and right-hand 
bubbles, respectively, at time t = to- To coincide 
with Eqs. (|7I8[) . the sign of has been chosen 

opposite to that of w^ =2 in one of the bubbles, while 
u>" =1 and w° =2 have been chosen to have the same 
sign in the other. The reason for the choice of the 
initial time derivatives is explained below. 

As discussed at the end of Sect. IIVE 11 at t = to 
and the functions wl{x) and w R (x) have the same 
magnitudes and shapes. Thus we write 



w L (x) = n w f w (t ,r,z + zo) 
w R (x) = n w f w (t ,r,z - zq) , 



(60) 



where f w (to, r, z ± zq) chosen similar to f s (to, r, z± 
Zo) of the scalar field so that the W fields fill the 
bubbles uniformly [12| ]. The normalization constant 
nw is determined as in Refs. [H],[l4| allowing in 
Eq. QUI) to be identified with w^ d rather than with 
. Of course for times t > to the distribution 
of the W in the bubble will differ from that of the 
scalar field since these fields evolve through different 
EOM. 

Boundary conditions on all the fields wj and their 
time derivatives are also required to solve the EOM. 



According to the requirements of Sect. lIVEl w a and 
Wo should be empty at t = to, so the boundary con- 
ditions on these quantities are 



w (ta,x) = w II (to,x) = 
Wo(t ,x) = wtfitoix) = . 



(61) 



Boundary conditions on their time derivatives are 
determined by the constrained initial conditions, dis- 
cussed below. 

The initial time derivatives in Eqs. (|58I59[) assure 
that W of the plasma displaced by the expanding 
bubble end up populating the coherent mode inside 
the bubble in our 0(3) formulation. To achieve this, 
we take 



dw L (x,y,z) 
dt 

dw R {r,z) 
dt 



i df w (t ,x,y, z + zq) 
lw dt 
, df w (t ,x,y,z - Zq) 

v dt 



(62) 



w is discussed in Appendix [B] 
n' w is determined by f w and 
. In general, n' w ^ nw, 



where the choice of n 
As long as v waU > v w 

f s according to Eq. ([j_. U |i . -m gcuciai, ,* w 
and the simple and natural choice f w = f s at t = t 
gives nyy = n' w . Perhaps unexpectedly, with the 
motion of the surface taken into account in our 0(3) 
formulation, the analog of the boundary condition 
on the derivative of the w z field in Eqs. (|7|8|) of the 
0(1, 2) formulation that maintains the average den- 



sity of the W is not n\ 



w 



0. 



With the motion of the surface taken into ac- 
count in this fashion, the boundary condition on 
the derivative of the field continues to be the mech- 
anism by which the average density of the W re- 
mains constant inside the bubble and expands with 
it. This happens in the 0(1,2) formulation of 
Ref. [12] as a consequence of the boundary condi- 
tion dw z (T , z)/dr = 0. 

Choosing the distribution in Eq. (1501) to have the 
shape of the scalar field f s (t,r,z) in Fig [1] the 
boundary conditions w£(to, r, z) and w 1 / (to, r, z) are 
shown in Figs. [5] and [3] (with nw = 1). We will see 
explicitly how the em currents and associated mag- 
netic fields begin to form in the collision once the 



Note that the speed of a W boson in the surface of the 
bubble is given by v w ks p w / J p^y + where pw ~ 

h/Ro s (to) is determined by the zero point motion of the 
W in the bubble and is the W mass in the surface of 
bubble. For cases of interest v-w is expected to be less than 
Vwall- 
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FIG. 2: The initial condition wi(to,r, z) for r = 0. 
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FIG. 3: The initial condition wi 1 {to, r, z) for r = 0. 



bubbles begin to overlap when we solve the EOM in 
Sect. El 

For these calculations, we find it convenient to 
normalize the i»° fields to unity inside the bubbles 
rather than to % and n' w as in Eqs. (|60I62[) . The 
normalization re-enters the calculation as a factor 
in the normalization of the em current as found in 
Scct. UVFl 



Auxiliary condition and constrained initial 
conditions 



F. Cylindrical coordinates 

In view of the axial symmetry of the collision as 
specified in Sect. IIVBI with cylindrical coordinates 
the natural coordinate system for expressing results, 
the EOM in Eq. (23 with w%(x) given in Eq. 03]) 
and em currents in Eq. (l30|) with j v (x) given in 
Eq. (1441) are expressed in cylindrical coordinates in 
Ref |14l |. For convenience, they are reproduced here, 
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di 



dz 



These EOM are equivalent to Eqs. (|20|) when 
constraints imposed by the auxiliary condition of 
Eq. (l22j) are satisfied; see Appendix [Clfor the details. 
When these constraints are satisfied, Eqs. (|63l64l65p 
determine the W fields, w". 

Equation (|4"4)) expressed in terms of these wj pro- 
vide the em currents given in Ref. [l4| , 



Ge ab3 (w b d-w a +w«^ 



+ r"w" 



dw b 

~dT 



, dw 



zb 



dt 



(66) 



We have used the freedom in choosing the bound- 
ary conditions to match our calculation to that in 
Refs. [H, However, as discussed earlier, the ap- 
pearance of a spatial-dependent mass m(x) means 
that not all of the initial conditions are independent 
and must satisfy the constraints in Eqs. (|27|28[) . Al- 
though this is a complication that was not a sig- 
nificant issue in Ref. (l2T |. in 0(3) symmetry they 
are easily maintained in cylindrical geometry, which 
is used for obtaining numerical results in subsequent 
sections. We will use these conditions to fix the time 
derivatives of w a and Wq at t = to using results found 
in Appendix [U] 



4irj(x) = Ge aM (w b d-w a +rw a 



dw b 
dr 



r dr 
„ a dw b 

2 <lft 







1 n dw zl 
r dr 



, dw 



^r- 2 <^r)> (67) 



dz 



4irj z (x) = Ge abA {w b d-w a 

„ dWn n „ 

- <i£ +r w ir z 



dz 



2wZ 



dw b 

~~dT 



2rw a 



dw b 
dw b 



dr 
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With these expressions, current conservation 



where 



d v j v {t,?,z) 



djo(t,f,z) 
dt 

dj z (t,f,z) 
dz 



1 dr 2 j(t, r, z) 
r 



dt 



= o 



(69) 



can be shown to be satisfied as expected, and 
Maxwell's equation for B, Eq. (156"]) . gives the de- 
sired seed fields. From these results we can infer 
the importance of surface dynamics for generating 
magnetic fields in bubble collisions. 

Adopting the convention that w% with BCI and 
BCII are normalized to unity (in units of raw) in 
the bubble at t = to, the overall normalization 
of the em current is the product of the factor G 
and the square of nw introduced and discussed in 
Sect. lIVE2l Thus 



G — > n w G , 



(70) 



where G is calculated using the values of g and g 1 
quoted below Eq. (J3J and the value of % is spec- 
ified in Eq. (72) of the Appendix of Ref. [3, fix- 
ing the average number density of W in the bubble 
equal to the number density of those W quanta in 
the thermal plasma that can make a transition into 
the bubble without violating energy conservation. In 
this way, we find that Eq. (j?0"f becomes 



' W 



G w -38.5 + 1.36T C GeV 

= 2.32m w , (71) 



taking the transition temperature (in GeV) to be 
f c = 166 from Ref. [13 ■ 

With the em current having the form given in 
Eq. (l4"4l , one can show that the solution of Eq. (f31)| 
takes the form 



B = 



r 



(72) 



where (— y,x, 0)/r is a unit vector in the azimuthal 
direction. Thus, the magnetic field generated in the 
bubble collision lies in the azimuthal plane, with 
x and y components only, and encircles the z-axis, 
just as in the Abelian Higgs model and our earlier 
work pjfl . 

Using Eqs. (|36IT2[) and including the conductivity 
current from Eq. (1321) . we immediately obtain B^ in 
cylindrical coordinates as the solution of 



, d 2 d 2 1 n 8 1 d 2 d . „ 



'dt 2 
4tt/ 



dt' 



(73) 



or dz 
and we have used the fact that 

dB 



d x E = - 



dt 



(74) 



(75) 



In obtaining the above results, the following iden- 
tities are useful, 

d-w a = -^- + 2w a + r— ^ , (76) 



dt 



dr dz 



wa - 9 = <l +rwa l- wa 4 z > (77) 



and 



w a ■ dm 2 Wq dm 2 rw a dm 2 w a z dm 2 
m 2 m 2 dt tit 2 dr m 2 dz 

„ , „ d In m n dlnm 
= 2«— — +rw a - 



dt 
, 9 In to 



dr 



- w, 



dz 



(78) 



V. NUMERICAL RESULTS, a s 



Because the magnetic field is critical for determin- 
ing whether present day galactic fields could have 
been seeded during the primordial EWPT, it is the 
most important prediction of our theory and the fo- 
cus of our numerical study. In the process of de- 
termining it we will explore the role played by the 
bubble surface in the production of these fields, as 
this has not been examined in previous studies. To 
achieve this understanding we will examine the W 
fields, the source of the em currents in our MSSM 
formulation, as well as the currents themselves. 

In this section we assume v wa u = 1. Although 
our theory as formulated accommodates wall speeds 
Vwaii < 1 and a ^ 0, we postpone discussion of these 
effects to a later section. 

Our scalar field in an isolated bubble for a s — 4 
is the same as the one in our example of Sect. IIV CI 
With the scalar field of an isolated bubble chosen 
in this way, the scalar field describing the collision is 
the same as the one shown in Fig.[TJ The collision as- 
sumed to take place in an average scalar background 
field of p av = -lmw (Pav/po — 0.046), correspond- 
ing to a collision early in the evolution of the phase 
transition. 
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With this scalar field, at t = to the radius of each 
bubble is -Ri/2^0) ~ 24, so the bubble surfaces are 
separated on the z axis by 2(z — -^1/2(^0)) ~ 22. 
Although we show results for to < t < t c + 20, 
we have solved the EOM over the longer interval 
to < t < t c + St max , where with 5t max = 29, at 
which time the radius of each bubble has increased 
to i? 1 / 2 (40) sa 64. The numerical accuracy of the 
solution deteriorates rapidly for 5t > 20. 

We show results for both nucleation and colli- 
sions. The nucleation stage corresponds to times 
t < t c = 11, during which time the bubbles may be 
considered as evolving approximately independent 
of one another. Collisions then correspond to the 
interval t c < t < t c + 20. 



A. Nucleation stage of bubble evolution, a s — 4 

During the nucleation stage of evolution, prior to 
the collision, both colliding bubbles have scalar fields 
with approximately the same shape and wall speed. 
In this case, the W a fields are proportional for a = I 
and a = II, and it is therefore sufficient to examine 
the fields in just one of the two bubbles. For this 
reason, for nucleation we drop the superscript dis- 
tinguishing the two boundary conditions. To sim- 
plify the discussion of nucleation, we calculate the 
W fields in a coordinate system translated so that 
nucleation originates at the origin, z = 0. 

The boundary conditions for our numerical simu- 
lations of nucleation are 

w (t ,x) = 
w(to,x) = 
w z (t ,r,z) = f w (t ,r,z) 

-w z (t ,r,z) = , (79) 

where f w (to,r, z) is the same quantity appearing in 
Eq. (pjUf . Because these are the bubbles that initiate 
the collisions presented in Sect. IV Bl below, the pa- 
rameters of f w (to j r, z) are clearly the same as those 
for the collision (with nyy = n' w ). Using the bound- 
ary conditions of Eq. (1791 . Eq. (|C17[) specifies the 
initial condition for dw(to,x)/dt and Eq. (|C18[) the 
initial condition for dwo(t , x)/dt. The W fields are 
then calculated solving the EOM given in Sect. lIVFl 
with boundary conditions specified in Eq. (|79[) with 
the scalar field of the bubble appearing in the exam- 
ple of Sect. IIV CI The resulting W fields are shown 
as a function of (t, z) for r w in Figs. I4I5I6I 




FIG. 4: Evolution of w z field in a bubble after nucleation. 
The field is plotted as a function of (t, z) for r ~ 0. 




FIG. 5: Evolution of w field in a bubble after nucleation. 
The field is plotted as a function of (t, z) for r 0. 



It is clear from the results that the W fields un- 
dergo a time-dependent evolution from nucleation at 
t = t„ = to the time of collision at t = t c = 11. 
The dominant component of the field is of course 
w z (t,r, z), since wq and w vanish at nucleation by 
choice of the boundary conditions. These smaller 
components of the field grow with time in part be- 
cause their time derivatives, given by the initial con- 
ditions, are finite at t = to. There is no magnetic 
field associated with a single bubble during this in- 
terval for reasons discussed in Sect. IIV E 11 However, 
there is a small field generated prior to t = t c in the 
collision of two bubbles as their surfaces, which have 
a finite diffuseness, begin to interpenetrate. 
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FIG. 6: Evolution of wo field in an isolated bubble. The 
field is plotted as a function of (i, z) for r ~ 0. 



B. Bubble collisions, a s = 4 

We turn our attention to collisions, the source 
of magnetic fields. Using the boundary conditions 
given in Eqs. (|58I59[) and shown in Fig. [2] and 
Fig. [3] and the boundary conditions in Eqs. (|6"Tj) . 
Eqs. (|C13IC14p give the constrained initial condi- 
tions for dw 11 (t ,x)/dt and dw 1 (to , x) / dt . Like- 
wise, initial conditions for Owq 1 (to , x) / dt and 
dwl(to,x)/dt are given in Eqs. (jC15IC16[) , 

With the scalar field describing the collision as 
shown in Fig. [TJ the profile of the region of bubble 
overlap in the collision at t — t c + 20 ~ 31 is shown 
in Fig. [7] At this time, the collision region extends 
along the z-axis from — i? 1 / 2 (31) — z — —91 < z < 
fl 1/2 (31) + *o = 91. 

The fields w%(t, r, z), Wo(t, r, z), and w a (t, r, z) for 
the collision are determined by solving the EOM 
given in Sect. lIVFl on the interval to < t < t c +St max , 
where to = and 5t rnax = 29. The results are illus- 
trated showing w\ in Fig. [5] plotted in the x — z plane 
at t = t c + St with St = 20. The region over which 
w z 7^ delineates the region i = c of bubble coa- 
lescence, with the bubble collision region clearly in 
evidence. 

The em current j z (t,r,z), j(t,r,z) and jo(t,r,z) 
given in Eqs. (|66l67l68p and calculated using these 
W fields at time t = t c + 20 is shown in Figs. 19110111"! 
respectively. The current is confined to a region that 
extends over a distance of ±40 into the x — z 
plane and a distance of ~ ±20 along the longitudinal 
direction. This region is quite comparable to the 
bubble overlap region at this time as shown in Fig.[7J 
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FIG. 7: Showing the projection of the surface S c of the 
collision region onto the x — z plane at time tj — t c + 20 
with v wa ii = 1. The surfaces S c is defined by Eq. (|38[) 
with Ry 2 (tf) = 55. 




FIG. 8: wi(t,x) in the x — z plane at t = t c + 20 for 
collisions of bubbles. 



In Fig. [12] we show the current j^(x) appearing in 
Eq. ([75} for the azimuthal magnetic field. 



C. Magnetic fields in bubble collisions, a s — 4 

The magnetic field calculated from Maxwell's 
Equation, Eq. (|73[) . with the em current shown above 
and with boundary conditions as discussed is shown 
in Fig. Q21 To facilitate comparison with the analo- 
gous results of Ref. [l2| with 0(1, 2) symmetry it has 
been plotted at intervals St following the onset of the 
collision comparable to the results given there. As 
expected, the magnetic field moves away from r = 
and increases in magnitude with t. 
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FIG. 9: Current 4nj z (t,x) in the x — z plane at t — 
t c + 20. 




FIG. 10: Current 4vrrj(t, x) in the x — z plane at t = 
t c + 20. 

The corresponding fields calculated from Ref. [l2j 
are shown in Fig. [14] Comparing Figs. I13|14l our 
magnetic field is about twice as large. As in the 
0(1,2) model, it is confined predominately to the 
region of bubble overlap shown in Fig. [7] The mag- 
netic field is also apparently smoothed out by the 
surface so that our result in Fig. [T3] does not show 
the oscillations apparent in Fig. [14] and fills the re- 
gion of bubble overlap more uniformly. 

Thus, bubble surface dynamics seems to produce 
fields significantly larger in scale as well as magni- 
tude. It is possible that bubble walls of even smaller 
surface thickness might grow even larger. We quan- 
tify this in Sect. IVII Calculations for the magnetic 
field including conductivity and wall speed v wa ii < 1 




FIG. 11: Current 4njo(t,x) in the x — z plane at t = 
t c + 20. 




FIG. 12: Current inj^(t,x) in the x — z plane at t = 
t c + 20. 

are given later in Sect. IVII1 

We close this section by illustrating the impor- 
tance of the boundary condition dw z (to)/dt. Re- 
sults for our present theory with n' w = are shown 
in Fig. Qjj] Comparing Figs. Q2] and [T5] it is seen 
that the magnitude and spatial scale of the mag- 
netic field both decrease by setting n' w = 0. This 
can be understood as follows. There are two points 
to be made. 

First, as noted in Sect. IIV E 21 in the absence of a 
mechanism for W of the plasma to enter the bubble 
as it expands (nw = 0), the initial speed of expan- 
sion vw(t = to) of the W field is determined by the 
zero point motion of the W in the bubble and its 
mass in the surface. With our choice of parameters, 
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FIG. 13: Magnitude of the azimuthal magnetic field in 
the transverse plane at z=0 as a function of r for a series 
of times t = t c + St, where St =5 (solid curve), 10 (short 
dash curve), 15 (medium dash curve) and, 20 (long dash 
curve). Field is calculated as in this work for a s — 4, 
n' w = 1, v wa a — 1, and a = 0. 
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FIG. 14: Magnetic field calculated in Ref. [12| in the 
transverse plane at z=0 as a function of r for times t = 
t c + St. Legend is the same as in Fig. 1131 



the momentum arising from the zero point motion 
is 



1 



ft 



Os 



1/24 w 0.042 



(80) 



and the value of its mass in the surface is mw = 
gPav/V2 ~ .0457. This suggests that as the expan- 
sion begins, 



v w (t = to) 



Pw 



0.68 



m 



w 



(81) 



which is both smaller than the speed of the bubble 
wall, v wa a = c, and close to the speed of expansion 



FIG. 15: Magnitude of the azimuthal magnetic field in 
the transverse plane at z=0 as a function of r. Legend 
is the same as in Fig. 1131 Field is calculated as in this 
work for a s — 4 and n' w = 0, v wa ii = 1, and a — 0. 



the magnetic field seen in Fig. [15] As time proceeds 
the spatial scale of the W field, and hence the em 
current, increasingly lags the growth of the bubble 
wall. Eventually, mw approaches mw and the W 
becomes even more confined to the interior of the 
bubble. The second point is that by choosing nw = 
1 so that the number of W populating the bubble 
grows in proportion to the volume displaced by the 
bubble, the em current naturally grows as well and 
leads to the larger magnitude seen in Fig. Q2] 

This discussion gives the rationale for the specu- 
lation made in Ref. [14j that the scale and magni- 
tude of the magnetic field might be sensitive both to 
the boundary condition for dw z /dt as well as to the 
steepness of the bubble surface. 



VI. SENSITIVITY OF MAGNETIC FIELDS 
TO BUBBLE WALL THICKNESS 

Since our theory accounts for surface dynamics, 
we expect its consequences to be quite different from 
earlier studies that did not address this aspect of the 
collisions. This expectation is based on the obser- 
vation [3] that when the surface is considered new 
terms appear in the EOM that manifest a strong sen- 
sitivity to the steepness of the scalar field in the bub- 
ble surface. Based on the observations given here, 
one might expect the importance of the surface to 
grow as the surface becomes steeper and, conversely, 
less striking for a more diffuse bubble surface. 

A quantitative measure of the sensitivity to the 
steepness of the scalar field in the surface may be ob- 
tained by comparing the calculation shown in Fig. 1131 
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to calculations with sharper walls. This sensitivity 
has, to our knowledge, never been studied quantita- 
tively in previous work. To obtain this quantitative 
measure, we compare here calculations with differ- 
ent surface diffuseness, a s = 4 and a s = 3. We show 
how the magnetic seed fields behave for a s — 2 by 
comparing to the results of Ref. 14j . 
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FIG. 16: Magnitude of the azimuthal magnetic field in 
the transverse plane at z=0 as a function of r for a series 
of times t = t c + St. Legend is the same as in Fig. 1131 
Field is calculated as in this work for a s — 3 and n' w — 1, 
Vwaii = 1, and a = 0. 

The calculation in Fig.[T|)]is identical to the one in 
Fig. EH but with a s — 3. Comparing the magnetic 
fields in these two figures confirms that the peak 
size of the seed field B^ Iax (t = t c + 20, a s )/m'y V in- 
creases as the wall becomes sharper. In particular, 
decreasing a s from a s = 4 to a s = 3 results in nearly 
doubling this peak field. 

What about the peak field for a s = 21 Comparing 
Fig. El (a s = 4) to Fig. 3 of Ref. 14] (a s = 2), 
both of which are calculated with nw — 0, we see 
that decreasing a s from a s = 4 to a s = 2 results 
in nearly a 5-fold increase in the peak field. Since 
B 



c (t = t c + 20, a s = 4)/m 2 w = 8 from Fig[H we 
conclude that with nw = 1 the peak field would be 



B 



c (t = t c + 20,a s = 2)/m 



40! 



VII. 



SENSITIVITY OF MAGNETIC FIELDS 
TO v waU AND CONDUCTIVITY 



Taking v wa u < 1 and a ^ would of course mod- 
ify the results found in the previous section. To 
quantify their effects in our model, we examine the 
case a s = 4, taking v wa ii — 1/2 [19], and the same 
boundary conditions, shown in Fig. 12131 




FIG. 17: Showing the projection of the surface S c of the 
collision region onto the x — z plane at time tj = t c + 20 
with v wa u — 1/2. The surfaces S c is defined by Eq. (|3Bp 
with Ri/ 2 (tf) = 55. 




FIG. 18: Magnitude of the azimuthal magnetic field in 
the transverse plane at z~0 as a function of r for a series 
of times t — t c + St. Legend is the same as in Fig. 1131 
Field is calculated as in this work for a s = 4 and n' w = 1, 
v W aii = 1/2, and a = 0. 



We first show the sensitivity to v wa u taking a = 0. 
Results are shown in Fig. [18] To interpret them, it 
is useful to refer to Fig. [TTl showing the region of 
bubble overlap at t = t c + 20 for v wa u = 1/2. There 
are several points to make. 

First, note that the scale of the magnetic field in 
Fig. [TH] exceeds that in Fig. Q21 This is because the 
magnetic field propagates with the speed of light and 
thus quickly escapes the bubble into the surrounding 
plasma. Secondly, we see that the peak magnetic 
field is larger than the peak magnetic field in Fig. (T3J 
The reason for this behavior is that the current is due 
to the W fields, which propagate in the plasma with 
an effective mass that is less than mw', hence, even 
though the bubble surface propagates with a speed 
less than the speed of light, the current outside the 
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FIG. 19: Magnitude of the azimuthal magnetic field in 
the transverse plane at z=0 as a function of r for a series 
of times t = t c + St. Legend is the same as in Fig. 1131 
Field is calculated as in this work for a s — 4 and n' w — 1, 
Vwaii = 1/2, and a = 6.7T. 



bubble actually has a greater speed, allowing a larger 
magnetic field there. 

Next we show in Fig. [19] the same calculation but 
with a = 6.7T as used in Ref @, taking T = T c = 
199GeV [13 ■ It is again useful to refer to Fig. Q3J 
as well as Fig. [T5] There are several points to note 
here as well. First, the scale of the magnetic field 
has decreased in comparison to that in Fig. [T5] This 
has occurred at the expense of the field outside the 
bubble, a consequence of the dissipation character- 
istic a 7^ 0. Thus, the field tends to be confined 
or "frozen" in the bubble. Both of these points are 
both consistent with the discussion of Sect. Ill CI 

The net effect of v wa u < 1 and a ^ is to pro- 
duce fields somewhat smaller than and of compara- 
ble spatial scale to those shown in Fig. (T3J which was 
calculated in the absence of these effects. The fact 
that there is such a small net effect is perhaps un- 
expected based on the experience with the Abelian 
Higgs model as discussed in Sect. Ill CI The larger 
size found here has been explained from the physics 
of the current arising from the charged W as opposed 
to that arising from the phase of the Higgs field: in 
the latter case, the current is strictly confined to the 
region within the bubbles, but in the former case the 
current is not so restricted. 



VIII. SUMMARY AND DISCUSSION 

Using EOM derived from the MSSM, we have ex- 
plored magnetic seed field creation during a first or- 
der primordial EWPT, focusing on the role that bub- 



ble surface dynamics plays in creating such fields by 
extending the results of Ref. [l4| ■ In our theory the 
charged gauge bosons W ± of the non- Abelian sec- 
tor and fermions, through the conductivity current, 
are the sources of the em current producing the seed 
fields. We develop a linearized version of the theory 
applicable for the case of gentle collisions, where the 
Higgs field is largely unperturbed by the collision 
process [121 ] . 

In the theory developed here, the W fields con- 
tributing to the em current evolve from initial con- 
ditions applied before the bubbles collide. These W 
fields initially fill the bubble uniformly and expand 
at the same speed as the scalar field of the bubble 
containing them, consistent with the physical pic- 
ture [12| adopted in earlier studies. 

Because of our particular focus on the dynamics 
of the bubble surface, the EOM must be solved nu- 
merically to obtain the W fields for nucleation and 
collisions. The magnetic field obtained from these 
solutions were found to be larger in both magni- 
tude and scale than the corresponding one embody- 
ing 0(1, 2) symmetry and jump boundary conditions 
at the time of collision [IJ, [l2|]. By comparing re- 
sults calculated with surfaces of different slopes we 
found the seed fields to be quite sensitive to the bub- 
ble surface, and we obtained a quantitative measure 
of the sensitivity to the steepness of the bubble sur- 
face. These results help understand why our results 
are larger in scale and magnitude than those ignor- 
ing surface dynamics. 

We have not attempted to determine the present 
day magnetic fields that are seeded by our fields 
generated during the EWPT. This is a complicated 
problem of plasma physics that has been studied ex- 
tensively. It is known that the characteristics of the 
primordial magnetic field are vastly modified during 
cosmic evolution to the present day. To model this 
evolution in magnetohydronamics, one first solves 
the equations for the period leading to the forma- 
tion of galaxies and galactic clusters considering all 
relevant dissipative processes such as viscosity, and 
then for the evolution of these structures including 
the possibility that they provide a large-scale dy- 
namo. Such studies have led to quantitative predic- 
tions for magnetic field energy and coherence length 
at the present epoch. The most recent of these [22j 
support the possibility that galactic cluster magnetic 
fields may be entirely primordial in origin. 

Our present results reinforce the hope that the 
EWPT is a promising source for production of seed 
fields in the MSSM for large-scale galactic and extra- 
galactic fields observed today. 
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Appendix A: The Scalar Field of a Bubble 

In this appendix we describe our scalar field p(x) 
for a single bubble. This field is is taken to have the 
form 

p(t,r,z) = p c f s (t,r,z) , (Al) 

where f s (t,r,z) gives its shape and time- 
dependence, p c = pit, 0, 0) is its value at the 
center of a bubble, and f s (t, 0,0) = 1 fixes its 
normalization. 

As discussed in Sect. Ill Al in the absence of a 
medium p(x) is an instanton solution of Coleman's 
equation, Eq. (|2"TT) . at nucleation t = t n . The mag- 
nitude of this solution is constant at p c — po to a 
high degree of accuracy throughout a region near 
the center of this bubble and then drops to zero at 
the surface. The scalar potential po is determined 
by mw from Eq. (fl7|) as 

m 2 w = \pl ■ (A2) 

Taking mw and g from Eq. we find from 

Eq. 422), 

Po « 176 GeV » 2.19m w . (A3) 

Within a medium, many-body corrections to the 
scalar field may be taken into account by introducing 
an effective scalar field. In the absence of a detailed 
understanding, f s (t, r, z) is taken to have the surface 
characteristics of the Coleman solution in the so- 
called thin-wall approximation, 

f (t,r,z) = -777V 1 - t aim ,(A4) 

J sit) a s 

where 

/ s m = l-tanh ~ J?0s ^ (A5) 
a s 



maintains the normalization and p c — > po — p av , 
where p av accounts for the presence of the other 
bubbles on the average as discussed in Sect. IIVDI 
Equation (IA4[) has the structure of the scalar field 
in Refs. [14[ and resembles the scalar field of Eq. (6) 
in Sect. V of Ref. [lOj. The functional form given in 
Eq. (|A4[) is used to characterize the distribution of 
other constituents of the bubble as well. 

In Eq. (|A4I) . i?(r, z) is the distance from the center 
of the bubble to any point x — (r, z) in the medium, 
with a s determining the fall off of the scalar field in 
the bubble surface. The function Ro s (t) specifies the 
time-dependence of the bubble radius, 

Ro s (t) = R t + v waU (t - t n ) , (A6) 

and is parametrized in terms of the asymptotic wall 
speed v wa ii and nucleation time t n of the bubble. 

Defining Rx/^it) as the radius at which the 
scalar field is half its central value, we find from 
Eq. (IA4IA5I) that 

i? 1/2 W = ^log(2 + exp^% (A7) 
2 a s 

The radius of the bubble at nucleation, r ns = 
Ri/2(t n ), is then 

r ns = ^log(2 + e X p^) , (A8) 

showing that Rt ~ r ns for large Rt/a s . 

With the nucleation centers for the bubbles lo- 
cated symmetrically on the z-axis at z = ±Zq, 
the collision time t c for bubbles nucleated simul- 
taneously may be found by solving the equation 
#1/2 (*c) = z Q . Taking R Qs {t) from Eq. (|A6]) . 

tc = t n + log[ —} . (A9) 

zv wa ii _2-)-e~ 

Appendix B: Boundary condition on derivative 

Of W z 

Requiring that the rate at which quanta enter 
the bubble be the same as the rate at which they 
leave the modes from which they originated in the 
thermal plasma, we have 

^.^bubblei^) = -Jj_N^ asma it) . (Bl) 

This condition maintains the average density of W 
inside the bubble roughly constant as a function of 
time. 
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The quantity N$ asma is given by N^ asma = 
Pw{T)V bubb ie, where Vbubble is the volume of the 
bubble and pw{T) is the density of the W ± in the 
plasma at temperature T that are able to make the 
transition into the bubble. Thus, the right-hand side 
of Eq. (IBip may be written 

-TiNplasrnait) = PW (T) — Vbubble , (B2) 

where d/dtVbubbie is the rate at which the bubble 
displaces the plasma volume. 

Since the volume of the bubble is the same as the 
volume occupied by its scalar field, we may write 



VbubbUt) = / f s (t,xydv 



(B3) 



Then, taking p w (T) = 2m w n 2 v from Eq. (72) of the 
Appendix of Ref. [12j |. with yJ2mw the relativistic 
normalization of the W wave function, we find 2 

2m w n 2 w f dr{ t S)2 dV(m) 



P (T) — Vbubble 

dt 



Likewise, N™ bble is determined as 



dt 



N, 



w 



bubble V 



(t) = 2m w / w z (t,x) 2 dV . (B5) 



Thus, the left-hand side of Eq. (|Blj) may written 
d 



dt 



^bubble (*) 



w z (t,x) dv . (B6) 



dt 



or 



-Tj.NZbblM = ^n w n' w m w 



f w (t,x) dfW ^ ] dV .(B7) 

To obtain Eq. (|B7|) have used Eq. ([SO)) and taken the 
boundary condition on the time derivative of the w z 
field at t « to to be determined by the same function 
that determines the shape of w z (to,x) in Eq. (|60|) . 
given in Eq. (|rJ2")l . 

The normalization n^- in Eq. (|F37|) may now be 
determined by equating Eqs. (IB4IB7|) at t = to, with 
the result 

Jdf s2 (t Q ,x)/dtdV 



n w 



Jdf™ 2 (t ,x)/dtdV 



(B8) 



Appendix C: Constrained initial conditions in 
cylindrical coordinates 

In this Appendix we express the constrained ini- 
tial conditions discussed in Sect. lIVE 3l in cylindrical 
coordinates using the specific form of the W field at 
t = to motivated and discussed in Sect. IIV E 21 

Referring to Eqs. (|26I2TI28[) the auxiliary condi- 
tion in cylindrical coordinates is enforced by express- 
ing x a ( x )> defined in Eq. ([2"3")l . and dx a {x)/dt, with 
the aid of Eqs. (I76I78|) . This gives 

n , s dwf, „ „ dw a duo? 
x a (x) = -tS- +2w a +r- 



dt 



„ , „ d In to 
2{w%—— + rW 

„ d In to . 



dr dz 
„ d In to 



9r 



(CI) 



and 



<9x a (a0 

at 



9 /^tUn „ „ dw a dw?. 

= 7T+ ^f + 2w a^) 

at at or 02; 

„ d , . „ a In to „ a In to 



at" u dt 
5 In to.. 



3r 



(C2) 



Equation (|C2[) may be equivalently expressed as 

d x a (x) 



dt 



a a ^ dw a dw z d 2 WQ 
dt dr dz dr 2 



1 au>o d 2 WQ 



r dr 



dz 2 



(C3) 



by using the EOM, Eq. (|63|). and canceling the term 
d 2 wljdt 2 . 

Using Eqs. (|C1IC3I) . the conditions that x a (t = 
to, 2) = and dx a (t = to,x)/dt = become, re- 
spectively, 







, dw?; „ dw a dw°: . 
' 1 2w a - 1 ~ 



v at 

„ . a In to 
+ 2(Wq — — h nr 



dr dz 
, a In to 



at 

, a In to 

dz 



dr 



(C4) 



and 



2 To avoid confusing the notation for the magnitude of 
w z (tQ,x) with other components of the W field, we call 
this magnitude here nw (called w in Ref. [l2| ). 



= 2- 



a w Q a 2 w Q a 2 < 



at drdt dzdt 
'Wq 1 a^o d 2 WQ 



dr 2 



r dr 



dz 2 



m 2 w a . (C5) 
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In accord with Sect. IIV A[ at t = to for both nu- 
cleation and collisions, the W quanta occupy w° 
with the other components Wq and w a empty. Im- 
posing the boundary conditions w a (to,x) = and 
WQ(t ,x) = 0, Eqs. (|C4IC5I) simplify and become 

dw$(t ,£) dw%(t ,x) 







at 

2w a z (t ,x) 



and 



dz 

ain77i(to, x) 

dz 

32,,,a i 



(C6) 



drdt 



8w a (t ,x) d 2 w a (t ,x) 

d 2 w a z (t ,x) 
dzdt 



(C7) 



With m(t, a?) a continuous function, Eqs. (|C7|C6[) as 
well as the auxiliary condition Eq. (|2"21 are satisfied. 



1. Collisions 

The colliding bubbles, which are nucleated simul- 
taneously at time f = t„ at points z = ±zq located 
symmetrically about the origin on the z-axis with 
nucleation radii r ns , are naturally assumed to be well 
separated and non-overlapping before they collide, 
as required by the discussion in Sect. IIVBI Thus, 
the boundary conditions and initial conditions on 
the W fields are established well before the collision, 
satisfying the requirements of Sect. IIV El 

The functions Wl, wr and their time deriva- 
tives depend on r and z entirely through R = 
\J r 2 + (z ± zq) 2 at t = to, just as for the scalar field 
(see Eqs. (joTM^ ). Thus, 

dw z ±(r,z) z ± zq d 



dz 



R OR 

Z ± Zq 3 



W z ± 



a 



-w z ± 



(C8) 



and 



d w z ±(r,z) z ± zq d dw z ± 
= R ~d~R dt 

Z ± Zq 3 8w z ± 



dz 



dt 



r dr dt ■ < C9 > 
where the positive sign corresponds to the left-hand 
bubble and the negative sign the right-hand bubble. 

The initial condition on dw a /dt is found from the 
auxiliary condition in Eq. (|C7|) . 



(2 + r 



a . dw a 



i a 9 dw a 



dr dt 



r dr dt 
a dui^ 



dz dt 



(CIO) 



which determines dw a (t — to,x)/dt in terms of 
dw z (t = to, x)/dt. Using Eq. (|C9|) and the boundary 
conditions on w^(t ,x) in Eqs. (|58I59|) . Eq. (IUT0|) 
becomes 

ld_ r2 duf_ = z ± zq a dw a z± 
r dr dt r dr dt 

This may be integrated immediately to give 



d_ 

at 



w a (t Q ,x) 



z±z ,d a 



d_ 

at 



< ± (r = 0,z)). (C12) 



For a —II we find 

d_ 

at 



iif. -t\ z + z , d 
w (t , a;) = — 5— (— w L ir, z) 



r 2 dt 



at 



w L (r = 0,z)) 



z - z , a 



r 2 dt 



d_ 

at 



t«ji(r = 0,z)) (C13) 



and for a =1 
3 I 



at 



w (t ,x) 



z + z a 
a 

z- zo , a 

- J^ fi (r = 0,z)) . (C14) 



Next, using Eq. (|C6|) and once again the bound- 
ary conditions on w z (to,x) in Eqs. (|58I59[) . we find 
dwo/dt for BCII to be 



a^o (to,x) a 



at 



and for BCI to be 



d: 



wl J (t ,x)+2wl\to,x) 



a 

x — lnm(io,x) (C15) 
dz 



dwn(to,x) a t, „ r . 

°L = 7 -w I z (t ,x) + 2w I z (to,x) 
dt dz 

a 

x — lnm(t , x) . (C16) 



2. Nucleation 

Since the positively and negatively charged W 
field evolve exactly the same in a single bubble with 
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our initial conditions, the distinction between the Finally, using Eq. (|C6I) and the boundary condi- 

charged fields is immaterial for nucleation and we tions on w z (to,x) in Eq. (179p . we find dwo/dt as 

drop the superscript on w a (t 7 r, z) that distinguishes in Eqs. (|C15IC16I) . 

between W + and W~ . Then, assuming that the 

bubble is nucleated at the origin, r — z = 0, with 

nucleation radius r ns , and assuming that the scalar 

field of the bubble expands in accord with Eq. (|A6|) . 

the initial conditions for t n ~ to « t c are found as _ , _,. _ 

follows. dw^x)_ = d itQtjg) + 2w , M 

For dw/dt, the initial condition is found from oz 
Eq. (|C7|) with the boundary conditions on w z (to,x) x — lnmftn x) (C18) 

in Eqs. $7S§. As in Eqs. (IC13IC14|) . we find dz 

z 

-gm{t ,g) = —(—w z (r,z) 

- ^w z (r = 0,z)) . (C17) 
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